de Haas-van Alphen effect in single crystal MgB 2 
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We report observations of quantum oscillations in single crystals of the high temperature super- 
conductor MgE>2. Three de Haas- van Alphen frequencies are clearly resolved. Comparison with 
band structure calculations strongly suggests that two of these come from a single warped Fermi 
surface tube along the c direction, and that the third arises from cylindrical sections of an in-plane 
honeycomb network. The measured values of the effective mass range from 0.44 — 0.68m e . By com- 
paring these with band masses calculated recently by three groups, we find that the electron-phonon 
coupling strength A, is a factor ~ 3 larger for the c-axis tube orbits than for the in-plane network 
orbit, in accord with recent microscopic calculations. 



The discovery of high temperature superconductivity 
in MgB 2 0| has led to great interest in this hexagonal lay- 
ered compound. Early isotope effect work clearly showed 
that phonons are important Q, but subsequent experi- 
ments and calculations by a number of groups have re- 
vealed some unusual features, such as the possible exis- 
tence of two distinct superconducting gaps || ||, 5j, |f3|, 
that may be associated with two different bands [kTH. 

The bandstructure of MgB2 has been calculated by a 
number of groups |9) and recently, photoemission experi- 
ments jlO| have confirmed its qualitative features. Stud- 
ies of the de Haas- van Alphen (dHvA) effect provide cru- 
cial information about the electronic properties of metals 
and superconductors. Observation of quantum oscilla- 
tions allows the shape of the Fermi surface and the effec- 
tive masses of carriers on individual Fermi sheets to be 
found. Then, together with band-structure calculations, 
electron-phonon coupling constants can be obtained [jll) . 
Here we report a detailed study of the dHvA effect in sin- 
gle crystals of MgP>2. Three dHvA frequencies are clearly 
resolved in our data, and can be assigned to two distinct, 
orthogonal Fermi surfaces parallel to the main symme- 
try axes of MgB2 ||. The effective masses correspond- 
ing to these three frequencies have been measured and 
compared with recent calculations by three groups. The 
comparison shows clearly that the electron phonon en- 
hancement is large for the c-axis tube and much smaller 
for the in-plane electron-like tube. 

The single crystals used in this work were grown in 
Tokyo by high pressure synthesis, as described in Ref. 
Jr2[ . Two crystals from the same batch were studied in 
parallel using sensitive piezo-resistive cantilevers to mea- 
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FIG. 1: Panels a-c: dHvA signals at several angles, rotating 
from [001] to [100]. Panel d: FFTs of the data at 9 = 20° (Fi 
and F2) and 70° (-P3). The top and bottom panels show data 
for crystals A and B at 1.35 K and 0.36 K, respectively. 



dimensions of the crystals A and B were 230 x 80 x 20/im 3 
and 230 x 200 x 47^m 3 respectively (the shortest side be- 
ing along the c-axis) . Thermal contact between the crys- 
tals and nearby calibrated thermometers was ensured by 
immersion in 3 He or 4 He liquids or a few mbars pres- 
sure of exchange gas. Both cryostats had a single-axis 
low-temperature rotatable sample stage and angles were 
measured to a relative accuracy of at least 1°. For crys- 
tal A, the alignment was verified from the symmetry of 
high resolution, torque versus angle sweeps near the o 
axis at 15 T, while for B a correction of 6° was applied to 
make the dHvA frequencies symmetrical about the crys- 
tal axes. The data obtained in both sets of experiments 
were essentially identical. 
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Fig. Jl (panels a-c) show torque versus magnetic field 
data p5 for both crystals as a function of angle, as the 
crystals were rotated from B || [001] (9 = 0°) to [100] 
(6 = 90°). Initially, the orientation of the c-axis was 
determined from the platelet shape and the anisotropy 
of the upper critical field H C 2- As expected, the torque 
signal goes to zero at the principal symmetry directions 
(6 = 0° and 90°) but is already sizable only a few degrees 
away. The signal was easily resolved over a wide field 
range, giving well defined dHvA frequencies [see the fast 
Fourier transforms (FFTs) in Fig. [I] (panel d)]. Two 
frequencies (Fi and F%) were observed for field sweeps 
carried out within ^45° of the c-axis [Fig. 0(a)], while 
only a single frequency F3 was found for fields nearer the 
plane [Fig. 0(b) and (c)]. 

The dHvA signals were analyzed in the conventional 
way using the Lifshitz-Kosevich (LK) expression (ll], |l^] 
for the oscillatory magnetization of a 3D Fermi liquid 
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where F is the dHvA frequency [F — (Ti/2ne)A, A is 
the extremal orbit area in k-space]; 7 is the phase; i?£>, 
Rt and Rs are the damping factors from impurity scat- 
tering, temperature (T) and spin splitting respectively. 
Rd = cxp(— mriB / eBr), where tub is the unenhanced or 
'bare' band mass [[ll], |l6| and r is the scattering time. 
R T = X/(smhX) where X = (m* / m e ) (14,69 /B)T, and 
to* is the quasi-particle effective mass, that is enhanced 
over to 5 by both electron-electron and electron-phonon 
interactions, and m e is the free-electron mass. The spin 
splitting phase factor is given by Rs = cos(mrgms /2m e ) 
where g is the Lande <?-factor for electrons in a given 
orbit, and ms/ms is the enhancement of the Pauli sus- 
ceptibility from electron-electron interactions alone. 

Fig. |^ (upper panel) shows the angular dependence of 
the three dHvA frequencies observed in both crystals on 
rotating from B\\ [001] to [100]. There is excellent agree- 
ment between the two sets of data. The strong, nearly co- 
sine angular dependence of F\ and F2 suggests that they 
arise from cylindrical sections along the c-axis. The small 
deviations from a strict 1 / cos 9 dependence are shown 
in inset (a) of Fig. |[ Assuming a simple cosine c-axis 
dispersion [i.e., Sk = fo 2 (k 2 + ky)/2m — 2t c cos(cfcz)], Ya- 
maji [fl7| derived a formula which accounts for the 'magic- 
angle' magnetoresistance maxima in quasi-2D organic su- 
perconductors. This gives the angular dependence of the 
two frequencies, F± , arising from extremal areas of a sin- 
gle warped cylinder as, 

F±(6) cos = (h/2ne)[Kk) ±4.izmt c Jo(ckf ta-a.9)} (2) 

In this expression, the first term is the mean frequency 
(h/2ne)TTk 2 f = (ff + F§)/2, the prefactor of the second 
term (T%/2-Ke)4TTmt c is (F® — F®)/2 and Jo is the Bcsscl 
function. Here Ff and F® represent F\ and F2 at 9 = 0, 
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FIG. 2: Upper Panel: Angular dependence of all three 
frequencies for a rotation between [001] and [100]. The 
open/closed symbols are data for crystals A/B respectively. 
Inset (a) F\ cos 6 and F2 cos 6, with a fit to Eq. (Q). Inset (b) 
F3 cos(# — 90). Lower Panel: In-plane rotation plot for F$, 
crystal A, with fits: 2685 / ( | cos /3 sin <f>\ + | sin f3 cos 15 cos <f>\) 
where (3 = 60 or 120. 



and may be taken directly from the data, so this equation 
has no free parameters. The good agreement shown in 
inset (a) of Fig. ^ implies that F\ and F 2 are the two 
extremal orbits of a single warped tube. As shown in 
inset (b), for the other orbit, F3 cos(90 — 9) is remarkably 
constant (to within ±1% over ± 40°) implying that it is 
cylindrical with very little warping. 

In Fig. H (lower panel), we show results obtained in a 
second experiment on crystal A in which the c-axis was 
aligned at ^15° from the rotation axis and the magnetic 
field was rotated approximately within the basal plane. 
Two minima separated by 60° were observed, consistent 
with the hexagonal symmetry of MgB2. The minimum 
value of F3 in the upper panel agrees with that in the 
lower panel (after correcting for the 15° offset) showing 
that when 9 — 90°, in the upper panel, B lies along the 
[100] symmetry direction. 

Temperature dependent studies were made for all three 
frequencies from 1.35 K to 12 K (crystal A) and from 
0.36 K to 12 K (crystal B) and analyzed in the standard 
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TABLE I: Summary of results for both samples along with 
estimates of the frequencies and band masses from band- 
structure calculations (c) [ ppf . 
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way by fitting the data to the damping factor Rt [Eq. 
(Q)]. Fits for crystal B are shown in Fig. |L and values 
of to* obtained in this way for both crystals are given in 
Table Measurements of to* were made with B approx- 
imately 20° off a symmetry axis. The values in Table | 
have been multiplied by cos(0) [or sin(0)] to correct for 
the usual angular dependence of the mass for a tubular 
band. This correction is less than 8%. 

In the inset to Fig. || we show the angular dependence 
of the dHvA amplitude near [100], for both crystals. 
There is a pronounced dip between 14° and 18°. The 
slight difference for the two crystals is ascribed to small 
alignment errors about an axis perpendicular to the axis 
of rotation. No such dip was observed for F% and F2. 
We believe the dip for F3 is a 'spin-zero', which is of- 
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FIG. 3: Dependence of dHvA amplitude on temperature. The 
amplitude was measured by fitting the sine function [Eq. (1)] 
to 1.5 oscillations in the range 17.70 — * 17.82 T for F3 and 
16.6 — > 17.9T for F\ and F2. The field was determined as 
the mean of 1/B in this window. The fit is to Rt in Eq. 
([!]). Inset: Dependence of amplitude on angle for 8 near [100] 
for crystals A (o) and B (•). Raw data are shown for A at 
15 T and 1.35 K Q with data for B at 18 T and 0.36 K, 
normalized to the same peak height. 
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FIG. 4: dHvA amplitude (A s ) at 0.36 K versus inverse field 
for several values of 8. The inset shows up and down field 
sweeps at 81.7°. The field sweep direction is indicated. The 
vertical arrow marks the point where the hysteresis vanishes. 



ten observed in dHvA studies when the areas of spin-up 
and spin-down extremal orbits differ by a half-integral 
number of Landau quanta, leading to destructive inter- 
ference in the oscillatory magnetization. This dip was 
shown to persist to at least 4.2 K. Since F3 comes from 
a tubular surface, whose effective area oc 1/ cos 6*, it is 
easily shown that such zeros occur at angles 9 SZ given 
by o— = ncosf?,,, with n = 1,3,5... As before, — 
includes band effects and electron-electron enhancement 
but not electron-phonon enhancement of the standard 
type This is because the Pauli susceptibility of a 

metal can be enhanced by electron-electron interactions, 
e.g., by the Stoner mechanism, but not by the electron- 
phonon interaction. Our data give 9 SZ = 18 ±3°, and for 
n = 1 and g = 2 this gives = 0.476 ± 0.007TO e . 

The scattering rates of the three orbits were extracted 
using two procedures. We either subtracted a smoothly 
varying background and fitted the whole field sweep to 
Eq. (0) (or more precisely, T cx BM osc d\nF/d8), or fit- 
ted Eq. (0) to 1.5 oscillations in 1/B to give the dHvA 
amplitude with the thermal damping factor removed 
(A s ), as a function of B. In both cases we used the val- 
ues of to* from Table | in the Rt factor and set Rg= 1- 
In field regions where the sample is clearly in the normal 
state, the first method gives excellent fits to Eq. (Q) for 
all three frequencies. Scattering times derived from fits 
for both crystals and all frequencies are shown in Table 
P. Making the approximation that all three frequencies 
arise from circular areas in k-space, we can obtain kp for 
each area and hence vf = fikp/m*, and the mean free 
path £ = vpr. These are also shown in Table |. 

The second procedure removes thermal damping ef- 
fects and reveals changes in the amplitude (A s ) versus 
1/B plots caused by the onset of superconductivity be- 
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low H C 2. As shown in Fig. [| at angles of 70° or less, the 
usual exponential behavior is obtained, (with slopes oc 
l/cos(90 — 8) as expected for a tubular surface). How- 
ever, for angles nearer to the plane, the Dingle plots be- 
come non-linear because of additional damping caused 
by the growth of the superconducting gap, as has been 
observed in many superconductors (see for example Ref. 
p8[). This observation proves that our dHvA signals 
arise from MgB2 rather than any impurity phase. The 
inset to Fig. [| shows hysteresis arising from supercon- 
ductivity that roughly corresponds to the onset of non- 
linearity in the Dingle plot. However, the slope of the 
Dingle plot for 9 = 81.7° remains high above 17 T, sug- 
gesting that the superconductivity has not been entirely 
suppressed in the reversible region. 

Recently, detailed bandstructure calculations, includ- 
ing estimates of dHvA frequencies and masses, have been 
carried out by three independent groups 19, 2fj, M. All 
calculations predict two warped cylinders along the c di- 
rection and two honeycomb networks in the basal plane. 
Comparison of the calculated frequencies with our exper- 
imental data shows that F\ and F 2 arise from the smaller 
warped cylinder along the c-axis while F3 corresponds to 
the electron-like, in-plane tubular network whose median 
plane contains the A,H and L symmetry points. The dis- 
crepancies with theory are less than 300 T (see Table §) 
which is only 0.2 % of the area of the hexagonal Bril- 
louin zone. Therefore, we can confidently compare our 
experimental values of m* with the calculated values of 
rriB to obtain the electron-phonon enhancement factors 
(A). We find A = {m*/m B - 1)= 1.20, 1.08 and 0.40 
for Fy, F2 and F3 respectively, which compare favorably 
with the calculated values A = 1.25,1.25 and 0.47 (Refs. 

The detailed angular dependence of Fi , F 2 and F 3 has 
been calculated by Harima @. We find this also agrees 
with our data and supports the Yamaji analysis given 
above. The calculated values of tub for F 3 also allow us 
to estimate the Stoner enhancement of the susceptibility. 
We find the enhancement ms/niB — 1 = 0.51 on this 



orbit, which is twice the calculated value of 0.26 |20 



The remaining discrepancy with theory concerns the 
absence of other dHvA orbits in the present study, partic- 
ularly those arising from the second c-axis tube and the 
other in-plane network (band 4 in the notation of Ref. 
p9|). In view of the discrepancies of 300 T mentioned 
above, it is possible that the smallest orbits (< 500 T) are 
actually even smaller or absent. The non-observation of 
the other orbits could well be due to the relatively short 
mean-free-path of our crystals. Experiments to higher 
field and/or with purer crystals will clarify this point. 

In summary, we have presented dHvA data for two 
single crystals of the new superconductor MgB2 that 
are in excellent agreement with each other and in good 
agreement with the most recent band structure calcula- 
tions. The present work provides direct evidence that the 



electron-phonon interaction is large on the inner c-axis 
cylinder and much smaller on the in-plane honeycomb 
network. This supports microscopic theories of super- 
conductivity that invoke two bands with very different 
properties [m, ||. 
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